Asymptotic function for multi-growth surfaces using power-law noise 
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Numerical simulations are used to investigate the multiamne exponent a q and multi-growth ex- 
ponent f3 q of ballistic deposition growth for noise obeying a power-law distribution. The simulated 
values of f3 q are compared with the asymptotic function f3 q = | that is approximated from the 
power-law behavior of the distribution of height differences over time. They are in good agreement 
for large q. The simulated a q is found in the range - < ct q < This implies that large rare 

events tend to break the KPZ universality scaling-law at higher order q. 
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I. INTRODUCTION 

Growing rough surfaces occur everywhere in nature 
and are encountered in engineering and everyday life. Ex- 
amples include, fluid displacement in porous media [Q, 
the growth of crystals Q or colonies of bacteria [gj , the 
propagation of a wet front on paper Q, and so on ||. 
A growing rough surface is one of the simplest patterns 
created by a nonequilibrium state. Consequently, many 
studies have examined growing rough surfaces. Two ex- 
ponents characterize a growing rough surface. One is 
the roughness exponent a, which relates the space length 
scale x to the surface width it), [w ~ x a ). The other is 
the growth exponent (3, which relates the time scale t to 
the surface width, (w ~ t@). A scaling function includ- 
ing these two exponents can be written as w ~ x a ^( -^-), 
where z = ^ is the dynamic exponent || . There are two 
different scaling regimes in ^ depending on the argument 
u = \. ^(u) ~ vP when u <C 1, and #(u) ~ const, when 
u » 1. 

The KPZ (Kardar-Parisi-Zhang) equation was pro- 
posed as an equation that models the dynamics of a 
growing rough surface [Q. The KPZ equation is writ- 
ten as d t h = aV 2 h + ^(\7h) 2 + r), where h, a, A, and 
77 are the surface height, effective surface tension, lateral 
growth strength, and a noise term, respectively. Since 
the KPZ equation is a stochastic nonlinear differential 
equation, an exact solution cannot be obtained. How- 
ever, the roughness and growth exponents can be cal- 
culated by the renormalization group method as a — | 
and /3 = i, respectively 0. It has been suggested that 
the scaling-law a + z = 2 holds in the KPZ universality 
class. Although most numerical simulations follow this 
result, many experiments of growing rough surfaces show 
a ~ 0.75 — 0.85 > | §. However, the KPZ universality 
scaling law is satisfied even in many experiments whose a 
is larger than |. In order to interpret the large a, several 
models have been proposed. These include the power-law 



noise model j§K quenched noise model |9j , and correlated 
noise model JlOj] . Moreover, certain systems with a large 
a obey EW (Edwards- Wilkinson) universality O]. In 
this paper, we focus on the power-law noise model. 

The noise in KPZ growth is ordinarily considered un- 
correlated Gaussian noise. Zhang suggested KPZ growth 
with uncorrelated power-law noise [^j . He performed nu- 
merical simulations of a BD (ballistic deposition) model 
and found that a varies with the exponent of the power- 
law noise fj,. Power-law behavior of the noise distribu- 
tion was observed in a fluid flow experiment |l2| . There- 
fore, the Zhang model is considered a moderate model 
for growing rough surfaces. 

On the other hand, Barabasi et al. investigated multi- 
affinity of the BD model with power-law noise [|l3) . Mul- 
tiaffine analysis is defined by the gth order height-height 
correlation function C q (x) as C q (x) = (\h(x') — h(x' + 
x)\ q ) x > ~ x qaq , where a q is the qth order roughness ex- 
ponent. 

While the entropy spectrum method describes the dy- 
namic characteristics of self-similar fractals Jlif , it de- 
scribes the static characteristics of self-affine fractals [ p~5| . 
That is, multiaffine analysis is equivalent to the entropy 
spectrum method for self-affine fractals Therefore, 
we must discuss the behavior of the qth order growth 
exponent (3 q in order to examine the dynamic character- 
istics of growing self-affine fractals. In ordinary variance 
analysis for the case q — 2, a mean field approximation 
has been applied to obtain the relations among a, /3, and 
/i |lq| , and the scaling-law a + z = 2 was confirmed. 

Myllys et al. investigated the slow combustion of paper 
and reported that a q and fi q at the combustion front 
varied with q Ji7| . Based on their discussion, we can 
expand the qth order height-height correlation function 



C q (x, t) = (\5h(x' : - Sh(x' + x,t' + t)\") x ,^ , (1) 

where Sh(x,t) = h(x,t) — (h(t)) x . Then, we can define 
a q and (3 q as 
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C q (x,0)~!t"*<, 



(2a) 



2 



C,(0,t)~t* , « 



(2b) 







The exponents a q and /3 g are denned at the limits x 
and t —> 0, respectively. 

The above-mentioned numerical simulations did not 
examine the q dependences of a q and f3 q on various /i 
thoroughly. Barabasi studied the temporal fluctuation 
of the surface width, however, his analysis examined the 
period after the saturation of surface growth [|l8| . There- 
fore, we investigate the behavior of a q and j3 q directly to 
determine how a q and (3 q depend on the noise parameter 
fi for higher order q. We are also interested in whether 
the KPZ universality scaling law a + z — 2 is still valid 
for higher order q. Here, we report the results of numeri- 
cal simulations of the BD model, and compare simulated 
data and the calculated asymptotic functions for a q and 

/V 



II. SIMULATIONS 

We investigate a 1 -I- 1 dimensional BD model whose 
growing dynamics are described by 

h(x,t + l) = max.[h(x — 1, t) + r]{x — 1, t), h(x,t) 

+r)(x,t),h(x + l,t) + r)(x + l,t)]. (3) 

This ultra-discrete BD algorithm can be connected to the 
KPZ equation Jl9[] . We start with /i(x,0) = for all x, 
and the surface evolves according to Eq. We use peri- 
odic boundary conditions for the space dimension. The 
uncorrelated noise r\ is taken from a power-law distribu- 
tion in the form [pi 



V 



(m+i) 



for rj > 1; P(f]) — otherwise. 



We focus on the range 2 < /i < 5. The variance of P(tj) 
is finite for /i > 2 and its statistical properties differ from 
those of Gaussian noise for /i < 5 || |l6| . 

First, we carry out a BD model simulation to inves- 
tigate the qth order height-height correlation functions. 
Figure |l| shows the behavior of C g (cc,0) and C g (0, t) at 
[i = 3.0. From bottom to top, the curves correspond 
to q = 1, 2, . . . , 10. In Fig. [l], we use the system size 
L = 1024, the time step T =(a)5000, (b)2000, and the 
average ensemble number N =(a)100, (b)5. The expo- 
nents a„ and P q are obtained from the slopes of the curves 
in Fig. fy. The values of a q and j3 q are defined as the limit 
of x — > and t — > 0, respectively (Eqs. ^a] and [S]). In 
addition, for regions of large t and x, conventional non- 
multiaffine scaling appears fl3|| . Therefore, the fitting 
regions are approximated as hit < 3 and In a; < 3. Fig. ^ 
shows the q dependences of a q and (3 q on various /i. For 
large q in Fig. ||(b), (3 q seems to collapse independently 
of fi. Since most of the error for each point is within the 
symbol used, there are no error bars in Fig. ^|. 
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FIG. 1: The gth order height-height correlation functions at 
M = 3.0. (a) i lnC,(a;,0) vs. lnx. (b) i In C q (0, t) vs. lni. 



q " i \~i -I \ i q 

In the figures, the curves correspond to q = 1, 2, . 
bottom to top. 
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When we use noise with a Gaussian distribution 
obtain the fully parallel curves in plot Fig. 0. This indi- 
cates that Gaussian noise and power-law noise affect the 
qth order moment analysis differently, and that the sys- 
tem size is sufficient to calculate the qth order moment. 

Next, we measure the amplitude of the effective noise 
rj m = Sh(x t + 1) — 6h(x,t) from the growth patterns 
directly fl7f| . Fig. || shows a log-log plot of the prob- 
ability distribution P(|7? m |) vs. |?7m|- The parameters 
used are L = 1024, T = 5000, and N = 100. P{\rj m \) 
clearly shows power-law behavior, but these exponents 
are slightly different from — (fj, + 1). Since these distri- 
butions obey the power-law, we denote the exponent as 
— v in this paper. This difference between 1) and 

—v might result from the lateral growth effect of the BD 
model, which corresponds to the nonlinear term in the 
KPZ equation. 
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FIG. 2: The q dependences of (a,)a q and (h)(3 q . The solid 
and broken lines correspond to the asymptotic functions - 

and — respectively. The inset in (b) is a log-log plot. 



III. ASYMPTOTIC FUNCTION 

Here, we calculate the multi-growth height-height cor- 
relation function C q (0,t) based on the multiaffine con- 
cept ^Oj. We discuss the growth path at some po- 
sition x': Sh(x',t). Fig. || shows a schematic image 
of the growth path. The growth path 5h(x',t) is the 
intersection with the 5h(x, t) surface at some position 
x = x'. We can generally normalize the time range 
and height to T and Sh max — 5h m i n , respectively, i.e., 
we consider the path 5h(x',t) for the range < t < 1 
and < Sh < 1. We can use the local growth expo- 
nent 7 to characterize the local singularity of the surface 
growth \Sh{x',t' +r) - Sh{x',t')\ ~ t 1 , where r = T -1 
and 7 > 0. The number of height difference segments of 
length I = \Sh(x' ,t' + t) — 5h(x' , t')\ in the growth path 
can be written as l~ u from the results of Fig. There- 



FIG. 3: The distribution of height differences between neigh- 
bour times. The relation between fi (input noise parameter) 
and v (measured noise parameter) is plotted in the inset. 



fore, the number of segments on growth path N^dj 
that have singularity exponents in the range (7,7 + dry) 
is found to scale with r as N T (j) ~ t~ 1v . We can obtain 
the height-height correlation function for the limit r — > 
as follows 



C q (0,r) 



1 

T 



(r^^ r ( 7 )p(7)d 7 



(4) 



Where, the function ^0(7) is a density function indepen- 
dent of t. For a continuous system, the integral in Eq. 
|] must be dominated by the value of 7 that minimizes 
1 + 7 <7 — 71/. Therefore, we replace 7 with the value 7(g) 
and compare the exponent with Eq. 2b, 



& = -+7(«)(l--)- 
q q' 



(5) 



Note that 7(9) decreases monotonically with increasing 
q, and 7(g) > In the limiting case g > 1, we 

assume that 7(5) vanishes faster than i. Finally, we 
obtain a simple asymptotic form of (3 q as 



0a = -, at q > 1. 

q 



(6) 



We plot this function in Fig. || as solid lines. While this 
function does not include any fitting parameter, the sim- 
ulated data agree with this function for large q in Fig. 
||(b). If we assume that the relation a q + |p = 2 holds 
for higher order q, the asymptotic form of a q can be cal- 
culated as: 



1' 



at q > 1. 



(7) 
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surface profiles 



growth path 5h(x',t) 



5h(x,t,) 




FIG. 4: A schematic image of the growth path. The surface 
profile Sh(x,t') corresponds to a snapshot of the surface at 
some time t' . The intersection between the rough surface 
Sh(x,t) and some position x — x' is the growth path. 



This function is plotted in Fig. ||(a) as a broken line. The 
simulated values of a q seem to distribute in the region 
\ < a q < ^ at large q. 



IV. DISCUSSION 

The contribution of large segments in the integral of 
Eq. |i] becomes dominant for higher order q. Since large 
segments are characterized by small 7 in our notation, 
7(g) decays rapidly due to the presence of large rare 
events with power-law distributed noise. Then, the ef- 
fect of fjh (or v) becomes negligible, as written in Eq. ||. 
Moreover, Eq. |5| becomes equivalent to Eq. at q v. 
Since the absolute value of the second term in Eq. || is not 
small, (3 q deviates from Eq. ^ at q < v. This qualitative 



Donds to the phase transi 
1 



The term - originates 



change around q ~ v corres 
tion point of Barabasi et al. 

solely from the normalization factor A in Eq. || The 
presence of large rare events is necessary, but the value 
of the exponent of noise distribution is not important for 
obtaining the asymptotic result Eq. [| 

In the early growth stage, we can confirm the KPZ 
universality scaling-law perfectly. For instance, we cal- 
culated a q at T = 2000 and found that all a q ap- 
proach the curve for Eq. [?] independent of /i. In addi- 
tion, a crossover to conventional non-multiaffine scaling 
at large x is clearly observed. We can also observe power- 
law-like tails of the distribution of the height difference 



5h(x + l,t) — Sh(x,t). Then, a q leads to the function | 
for the limit T — ■+ 00 using the same calculation as for (3 q . 
Since the non- uniformity of the power-law noise increases 
as fi decreases, the influence of large rare events domi- 
nates in the small fi system. Therefore, the smaller /1 
becomes, the nearer a q approaches |. Namely, a q obeys 
the KPZ universality scaling law in the early grow- 
ing stage, and obeys the rare event dominant behavior | 
in the fully developed stage. 

In the inset of Fig. ||(b), a slight discrepancy between 
the simulated data and the function | is observed. The 
discrepancy is due to the small, but finite, j(q) effect. 
Since the probability of finding a large event decreases as 
fi increases, the value 7(9) does not decay quickly when 
\i is large. Therefore, the statistical property for large /1 
deviates from the power-law noise case and approaches 
the Gaussian noise case . The range 2 < /i < 5 is con- 
sidered the power-law dominant range. In this sense, Eq. 
^ might be a limited approximation form. The precise 
correction to use for the range of small /i (in particular, 
the Levy distribution case \i < 2) remains unsolved. 

In a paper combustion front experiment, f3 q varied 
across the value 0.5, and a q did not fall beneath 0.5 [ p7[ . 
This tendency in f3 q is similar to our result. Myllys et 
al. found the effective power-law noise amplitude in their 
system to be 3.72 <v< 5.0. Our simulations include this 
range. Then, the behavior of (3 q inevitably approaches 
the form of Eq. |6[ However, the behavior of a q differs 
from our result. This means that the paper combustion 
front grows according to power-law noise, but breaks the 
KPZ universality in a different manner from our result. 
This experimental system might belong to another uni- 
versality class. In other experimental systems, qih order 
exponents should be measured in order to discuss univer- 
sality classification in detail. There have been no mul- 
tiaffine analyses of other models, such as the quenched 
noise, correlated noise, or EW class models. Both nu- 
merical and experimental studies are needed to further 
understand growing rough surfaces. 

In summary, we performed numerical simulations of 
the BD model with power-law distributed noise (2 < 
fx < 5). The simulated qth order growth exponent is in 
good agreement with the approximate asymptotic func- 
tion j3g = - for higher order q. Assuming the KPZ uni- 
versality scaling-law a q + ^f- = 2, we obtain a q = 
oi q is distributed in the region ^ < a q < j^j- This indi- 
cates that a q is determined by competition between large 
rare events and KPZ universality. 
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